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We begin with an overview of the experimental results for the temperature and doping depen- 
dences of the optical-conductivity spectral weight in cuprate superconductors across the whole phase 
diagram. Then we discuss recent attempts to explain the observed behavior of the spectral weight 
using reduced and full models with unconventional d x 2_ y 2 charge-density waves. 

PACS numbers: 71.10.-w, 74.25.Gz, 74.72.-h 

I. PRELIMINARIES 

The optical-conductivity sum rule plays an important role in the understanding of the physics of cuprate supercon- 
ductors. The behavior of the spectral weight in the pseudogap phase and when the system undergoes a superconducting 
transition has been the object of the experimental investigation in the last few years [1-8]. Considering that the anal- 
ysis of the experimental data and the extraction of the spectral weight involves a lot of theoretical assumptions, so 
that even opposite conclusions can be drawn based on the same data [5, 9, 10], we decided to describe the current 
experimental situation in the introductory Sec. II using the theoretical language from the very beginning, to avoid 
any possible confusion. We find that a clear perception of these assumptions is crucial for the correct interpretation 
of the experiments. 

In Sec. Ill we briefly summarize the recent history of unconventional density waves, and introduce the reduced 
d-density wave (DDW) model, which can be regarded as the charge density wave equivalent of the reduced BCS 
model. Then based on Rcfs. [11, 12] we discuss the optical conductivity and the sum rule for the reduced model, as 
an attempt to interpret some of the recent data on the sum rule of cuprates. In Sec. IV we discuss the limitations 
of the reduced model, as already pointed out clearly in Ref. [12], and consider the recent results obtained for a more 
generic microscopic Hamiltonian by Aristov and Zcyhcr [13]. Some of the remaining questions and the relationship 
to other systems are discussed in Sec. V. 

II. GENERAL NOTATIONS AND OVERVIEW OF EXPERIMENTS 

A. Optical conductivity and sum rule for it 

The optical-conductivity sum rule is intimately related to the general properties of the current-current response 
function [14, 15]. Since we are interested in describing lattice systems, let us consider a general case of the electronic 
system described by the Hamiltonian: 

H = H + H int , H = -^2,Ujc\ a c ja - fi^24 a c ir j, (2.1) 

ij i 

where the field operator ct creates an electron of spin a at the site r,, fjj is the hopping parameter, /i is the 
chemical potential and Hi nt is the interaction term. For the kinetic part Hq a proper description of the cuprates 
is obtained by restricting the sum to next-neighboring sites on a squared lattice that results in a band dispersion 
e(k) = — 2t(cosk x a + cos k y a) + At' cos k x a cos k y a, where a is the lattice constant. However, to avoid unnecessary 
complications we will also refer below to the simplest case t' = 0. Throughout the paper units ft = ks = c = 1 are 
chosen. 

Given the model (2.1), the definitions of the particle current j(rj) and the so-called diamagnetic tensor Tu(ri) follow 
from the first- and second-order derivatives of the Hamiltonian if (A) with respect to the external vector potential A 
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[16]: 



H(Ai) « (0) - J2 



eA i (T J )j i (T i )-—A$(T i )>ni(T J ) 



(2.2) 



so that the total electric current density is expressed as 

J;(r) = -5H/5Ai(r) = e i*( r ) - e 2 r il (r)A i (r), (2.3) 

By using the generalized notation j„ = (ji,jo) where the index /x = (i,0) with i = 1,2 indicates spatial and time 
components respectively, and jo is the particle density, one can evaluate ( J M (g)) within linear response theory [15-17] 
and obtain J^q) = e K fll/ (q)A I/ (q), where the electromagnetic kernel is defined as: 

Kpvfaitom) = -'WW 1 ~~ ^o) +n Mi/ (q,iO m ). (2.4) 
Here n M „(q, zf2 m ) is the current-current correlation function 

1 f 3 

U llv (q,iQ m ) = - J dTe^ T (r r ^(q,T)^(-q,0)), (2.5) 

t is the imaginary time, (3 = 1/T, Q m = 2irmT is the bosonic Matsubara frequency and N is the number of unit cells. 
Finally, using that A(ui) = E(uj)/i(uj + iO), where E is the electric field, one arrives at the famous Kubo formula 

. 2 Kn(q = 0,u)) . 2 < t u > -n„(q = 0,w) 

a w ) = -ie — — — — = le , (2.6) 

V(Lu + iO) V(u) + iO) 

where V is the unit-cell volume and the standard analytic continuation ifl m — > cu+iO was made. To avoid confusion, we 
will indicate the imaginary bosonic frequencies with ifl m , the imaginary fermionic frequencies with ioj n = in(2n + 1)T 
and the real frequencies with to. Since an isotropic system is considered we can omit the index i as done in the LHS 
of Eq. (2.6) and in what follows. 

To extract the real and imaginary parts of the optical conductivity defined in Eq. (2.6) we should consider the 
effect of disorder. Indeed, for the (ideal) case of a system without disorder by taking the real part of (2.6) one obtains 
Rea(uj) = (ne 2 /V)S(uj)[< t > — ReII(0,cj)] + (e 2 /V A )ImII(0, w)/u). The presence of a delta function, which means an 
ideal conductivity of the system, is an artifact of the assumption that no impurities are present. As soon as disorder 
is present one expects that in the normal state the following identity holds: 

Rcn(q = 0,uj^ 0) =< r >, (2.7) 

so that only the regular part of cr(w) survives, defined as usual as: 

_ , . e 2 Im±I(q=0,w) 

ReaM = 77 — • 2.8 

V (jj 

Thus, using the Kramers-Kronig (KK) relations for II(q = 0, u>) one can derive the well-know sum rule: 

W(T) = r Rea(^ = t- f 1 ^ = °^ dw = !^ Rc n( q = o,c = 0) = ^<r>. (2.9) 

Observe that Eq. (2.7) and (2.9) require that for the system in the presence of disorder, the dynamic (q = Q,u> — > 0) 
and static (oj = 0, q — > 0) limits of RcII(q, lu) commute. This is indeed the case for a disordered system, but not for 
a clean one, where usually Ren(0,cj) vanishes as ui — ► [16, 17]. 

As one can see from Eq. (2.9) the sum rule is directly defined by the diamagnetic tensor, which in turn depends, 
according to Eq. (2.2), on the way how the vector potential A enters the Hamiltonian of the system. When a 
continuum model is considered instead of Eq. (2.1), the kinetic term is expressed as / (— V) 2 /2m and A is inserted 
using the minimal coupling prescription — iV — > — iV — eA. For lattice systems the equivalent of the minimal coupling 
prescription is the so-called Peierls ansatz [14, 16, 18], which corresponds to inserting the gauge field A in Eq. (2.1) by 
means of the substitution Cj — * Cje _ie J A dr . In this case, it is clear that when the interaction term of the Hamiltonian 
is a density-density interaction, as in Eq. (3.2) below, only the kinetic hopping term is modified, while the interaction 
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term is gauge invariant (GI). As a result, the current/density operator and the diamagnctic tensor can be expressed 
(for small q) as: 

= ^H^( k ) C Lq/2 CT C k+q/2-' ( 2 ' 10 ) 
k,cr 

fo> = ^Ew^'' (2.11) 

k,(T 



where 



^(k) = «,l), (2.12) 
and (v£)i = dek/dki is the Fermi velocity. As a consequence, we obtain the restricted optical-conductivity sum rule: 

W(u M ,T) = r Rea u (co,T)dco ee W(T) = ^f(r H ) = ^ £ ^n k , a = (2.13) 

where the last equality, relating the spectral weight to the average value of the kinetic energy K, only holds for a 
nearest- neighbor tight-binding model (i' = 0), and n ( i is the dimensionality of the system. Observe also that here 
K includes the sum over two spin components, so that the factor 1/rid comes from the fact that J2k® 2£ k/dk 2 = 
2ta 2 cos kid = —(a 2 /rid) Sk £ k- ^ n writing Eq. (2.13) we introduced a cut-off wm in the frequency integration to 
make contact with the experiments, where only a limited range of frequency is accessible. However, the equivalence we 
established between Eq. (2.9) and Eq. (2.13) relies on the fact that even though in deriving Eq. (2.9) an integration 
up to an infinite cut-off is formally required, an intrinsic cut-off is provided by the energy u>m below which the 
tight-binding description (2.1) is valid. As we shall see, the definition of such a cut-off in cuprates is a quite delicate 
issue. Indeed, it is clear from Eq. (2.13) that the spectral weight depends both on temperature and interaction trough 
the occupation number 7ik,<r, and as a consequence it is a powerful tool to test theoretical predictions coming from 
different models. This behavior should be contrasted to the case of an electronic system described by a quadratic 
band dispersion eu = k 2 /2m. Indeed, in this case the tensor (r^) reduces to nSij/m, where n is the total carrier 
density, so that Eq. (2.9) reduces to the so-called f-sum rule 

o 

7rnp 

Rca{uj)duj = — — , (2.14) 

, m 

which is temperature and interaction independent. To make contact between the restricted and full f-sum rule it can 
be useful to define the frequency-dependent spectral weight: 

W(w,T)= Rc a(Lu',T)duj'. (2.15) 



In a lattice system one would expect that below lj = u>m only transitions within the conduction band are included, 
while at higher energies interband processes play also a role, 

u{lo) = a intra (u) + a lnt er(u), (2-16) 

such that <Jintra(u) = at lu > u>m and (Tinter(w) = at uj < u>m- As a consequence, the spectral weight W(ljm,T) 
only measure (7mtra(w) and gives a measure of the interactions acting on the system through the relation (2.13). When 
the integration frequency of Eq. (2.15) is extended above ujm one should evolve towards the f-sum rule (2.14), and 
the inclusion of the processes described by cr i „ ter (w) should give a temperature independent spectral weight. Observe 
that while many theoretical results have been provided in the literature about the behavior of the restricted spectral 
weight (2.13), based on different interacting model as Eq. (2.1), no clear theoretical understanding exists yet about 
the evolution of the spectral weight towards the f-sum rule as ui > com [14]. 



B. Conventional spectral- weight behavior 

The relation between the optical-conductivity spectral weight and the diamagnetic tensor (or kinetic energy) intro- 
duced in the previous section is particularly useful from the theoretical point of view, because it allows one to derive 
the restricted sum rule without doing an explicit computation of the optical conductivity. However, one can also be 
interested on the precise structure of cr(uj), and on the eventual transfer of spectral weight through a phase transition. 
Before discussing the experimental data in the cuprates it is worth showing some examples of the applicability of the 
previous formulas in simple paradigmatic cases. 
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• Temperature dependence in the non-interacting nearest-neighbor tight-binding system 

In this case, in Eq. (2.13) the occupation number is simply given by the Fermi function, = /(£k), where 
£k = £k — The main temperature dependence of the spectral weight (2.13) comes from the temperature 
smearing of the Fermi function, and can be easily evaluated using the Sommerfeld expansion. In = 2 
dimensions one obtains 

W{T) 1 ^ f W(0) 



(ire 2 a 2 /V) 



= -jfY, £ ^ = - / deN(s)ef(s fx) = j^j^ - c{»)T\ (2.17) 



where c(e) = (it 2 /6){eN'(s) + N(e)], and N(e) and N'(e) are the density of states (per spin) of the tight-binding 
dispersion and its derivative, respectively. For a flat density of states N = 1/2D where D is the semi- bandwidth, 
so that: 

TT 2 

W(T) = W{0) - —T 2 = W{0) - BT 2 . (2.18) 

Observe that for a layered system V — a 2 d, where d is the inter-plane distance, so that the quantity W — 
W/ (ire 2 a 2 /V) has the dimension of an energy (see discussion Sec. II-C below). The hopping t of the two- 
dimensional band dispersion can be related to the semi-bandwidth as D — At. Analogously, one can consider 
the tight-binding dispersion for t' = and estimate, at small electron density, N(e) = l/4wt. The slope 
c = 7r/24< so obtained is not far from the value of c obtained around half filling, see Ref. [12]. 

• Cut-off definition in the Drude model 

In metals displaying a conventional Fermi-liquid behavior the quasiparticlcs behave as free electrons with a 
renormalized mass m* and a characteristic scattering time r = 1/T. In the Drude approximation the real part 
of the optical conductivity reads: 

, ne 2 T/m* T W(T) 

1 + Ul' ! T Z TT 1 A + U! Z 

where we made the dependence on the lattice spectral weight W explicit. In Fig. 1 we report the integrated 
spectral weight as a function of the frequency. As one can see, at frequencies of the order of few times the inverse 
scattering time T, 90% of the total value W(T) is recovered, showing that in this case the cut-off frequency 
lum can be identified as ojm ~ 6r [4]. In a conventional metal one expects that Eq. (2.19) reproduces in good 
approximation the optical conductivity, eventually with a T slightly temperature dependent. Moreover, the 
Fermi-liquid paradigm ensures that the quasiparticlcs still have a Fermi surface, with an occupation number 
given by the Fermi function. The only effect of interactions is to renormalize the bandwidth, which can 
be for example determined with LDA calculations or extracted by some measurement of the Fermi surface 
(as for example by Angle Resolved Photoemission Spectroscopy (ARPES)). An example of such conventional 
Fermi- liquid behavior in gold has been indeed reported recently in Ref. [6]. 

To make contact with the notation that we will use below while discussing the experiments, we observe that for 
the Drude model (2.19) one usually expresses the spectral weight in terms of the energy Q 2 P = Airne 2 /to, which 
corresponds to the relations: 

<JDru d e = $^, W(U> M ) = ^. (2.20) 
I — ILO 4 

Observe that in the Drude model the scale fl p is also related to the zero of the real part of the dielectric function 
e, which is defined as: 

e(w =ei+«e 2 =£ 00 + (2.21) 

UJ 

where represents the screening by interband transitions. Using Eq. (2.20) we see that: 

e? rud ^)=eoo- f ^- 2 . (2.22) 

As a consequence, we see that when r £lp the real part of the dielectric function vanishes at a frequency loq 
given by: 

e? = - c 2 = J (l - JJ) , (2.23) 
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so that within the Drude model uq ps ftp/y^ can also be related to the restricted spectral weight. Moreover, 
since in conventional metals T <C Sip, in these materials the term "plasma edge" refers both to the edge of 
intra-band optical absorption and to the spectral weight itself of the intra-band conductivity. However, the 
temperature dependence of uq is directly related to the temperature dependence of Hp (or W) only if T is also 
temperature independent. Otherwise, a narrowing of the Drude peak with decreasing temperature can lead to 
a blue shift of uq which is not necessarily related to a change of the spectral weight W [5] . 



r=0.025D 

1 i 1 1 




' ' 1 1 ' 1 

0.1 0.2 0.3 0.4 0.5 

m/D 



FIG. 1: Integral of the optical conductivity of the Drude model N(io) = 2 Rc<70 rllt j e (Li/)da/, normalized to the total spectral 
weight W, as a function of uj/D, where D is the bandwidth. Inset (a); plot of the optical conductivity for W = 0.2D. Inset (b): 
real part of the dielectric function. Observe that here f!p = 2W = 0AD, where one can find the zero of the dielectric function. 



• Superconducting transition 

When a system displays a phase transition to a new ordered state, with CDW, SDW or SC formation, the 
occupation number changes temperature dependence with respect to the simple Fermi function characteristic 
of the normal state of Fermi-liquid systems. A typical example is provided by the SC transition, where within 
the BCS approximation one can show that the spectral weight (2.11) acquires the form: 



W(A,T) 



2ire 2 



E 



dkf 



1 



& tanh^ 
2T 



(2.24) 



where At = (Aq /2) (cos k x a — cos k y a) is the <i-wave SC gap and 



sc 



A? is the quasiparticle dispersion 



in the SC state. One can easily see that the spectral weight W(A, T) decreases as T is lowered below T c . Bearing 
in mind Eq. (2.13), this decrease can be understood as a consequence of the kinetic-energy increase in the SC 
state, where the states above fi become occupied due to the particle-hole mixing. 



C. The experimental data in cuprates 

From the experimental point of view, the determination of the cut-off % in Eq. (2.13) follows naturally from the 
observation of a minimum of the reflectivity (or of the conductivity) of cuprates at the so-called plasma edge up. 
This energy scale is of the order of ujp — 8000 — 10000 enr 1 in all the families of cuprates, and can be identified as 
the frequency uq at which the real part of the dielectric function crosses zero [6]. It is then usually assumed that 
transitions above up involve interband processes, so that up is a good candidate for the cut-off um used in Eq. (2.13) 
to relate the spectral weight W(um,T) to the diamagnctic tensor, according to Eq. (2.9). In this range one expects 
that the temperature dependence of W(T) gives information about the interactions acting on the system. 

While the identification between up and um is quite natural given the band structure of cuprates, the optical 
conductivity cr(u) does not display at all the conventional Drude-like behavior represented by Eq. (2.19). A detailed 
discussion of all the peculiarities of the optical conductivity of cuprates is beyond the scope of this paper, and can 
be found for example in the recent reviews [19-21]. We only recall here that quite often the optical conductivity 
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of cuprates is analyzed in terms of the so-called extended Drude model, which relies on an extension of Eq. (2.20) 
including a frequency-dependent optical mass enhancement m* — > m[l + A(cj)] and optical scattering time t(cj), 
according to the formula: 

a M = M ( f P 't+U TV ( 2 - 25 ) 
1/t(uj) — iu)[l + A(cj)J 

Although the representations of experimental data using Rccr(u;), Im<7(a;) and l/r(w) with 1 + X(uj) are formally 
equivalent, it has become rather popular to discuss the pseudogap behavior [19-21] in cuprate superconductors using 
the language of the optical scattering rate and mass enhancement. As we explained above, Qp is related the reduced 
spectral weight W. Indeed, normally Op is found from the sum rule (2.14) written in the form f£ M Recr(w) = Op/8. 
As stressed in [20] , the choice discussed here of the somehow arbitrary cut-off com introduces an uncertainty in the 
overall scale factor of the optical scattering rate and optical mass in Eq. (2.25). 

As one can see in Table I, Op has values around 2 eV in F^S^CaC^Os (Bi-2212), Bi2Sr2Ca2Cu30io (Bi-2223) and 
YBa2Cu306+ x (YBCO) compounds, and somewhat smaller values in La2- 2: Sr 2 ,Cu04 (LSCO) samples. The analysis 
of the optical data is usually performed by using a cut-off ojm temperature and doping independent, even though 
some variation in ojq has been reported for example in BSCCO compounds [1]. However, as stressed in Ref. [5], the 
temperature variations of lvq cannot be directly related to temperature variations of the spectral weight, because the 
temperature variations of uiq are also affected by the thermal correction induced by the scattering time, see Eq. (2.23). 
In Table I we summarize some values of the low-temperature spectral weight taken from the experiments of Refs. [1-8]. 
As we explained before, the strength of the restricted partial weight W can be expressed both as W = Op/4 or in 
terms of the kinetic energy K per Cu in-plane atom as in Eq. (2.13). To fix the energetic units, it can be useful to 
take into account the following relations. In most of the experiments the spectral weight W exp is defined as half of the 
value (2.13), because one integrates only the data at positive frequencies. Then from Eq. (2.13) in rid = 2 dimensions 
and (2.20) we have two relations: 

Op = AW = 8W exp = E C K, (2.26) 

where the energy scale E c in appropriate units is: 

27re 2 a 2 2ttc 2 27T_ „„ Tr * a ,„ 
E c = —— = — = —U.38eVA, 2.27 
V da 

and d is the distance between Cu02 planes. To compare the data between single-layer and multi-layer compounds 
we used for d the average distance between two neighboring layers, which corresponds approximately to 13.2 A, 6 A, 
7.5 A and 6.2 A in LSCO, YBCO, Bi-2212 and Bi-2223 compounds respectively. Since the spectral weight W exp is 
measured usually in _1 cm -2 one should also use the equivalence lcm _1 =0.21 fi _1 cm _1 [4]. Then one can obtain 
both the spectral weight in cm -2 (or equivalcntly in cV 2 ), and the kinetic energy K, given by 8W exp /E c = tip/E c . 
Finally, we note that the background dielectric constant Soo has typical values around 5 in cuprates, explaining the 
difference between ujp and Op in Table I. 



Ref. 


Compound 


ui M ~ oj P (eV) 


Op(T = 0) (eV) 


E c (eV) 


K{T = 0) (meV) 


[1] 


Bi-2212, T c =66 K, UD 
Bi-2212, T c =88 K, OVD 


1.2 


1.90 
2.04 


12.04 


300 
346 


[2, 3, 7] 


Bi-2212, T c =70 K, UD 
Bi-2212, T c =80 K, OPD 
Bi-2212, T c =63 K, OVD 


1 


2.02 
1.82 
2.12 


12.04 


339 
275 
373 


[8] 


Bi-2223, T c =110 K, OPD 


1 


2.0 


14.61 


276 


[4] 


YBCO, T c =57 K, UD 
YBCO, T c =91 K, OPD 


1 
1 


1.64 
2.07 


15.05 


180 

285 


[5] 


YBCO, T c =92.7 K, OPD 


0.5* 


2.04 


15.05 


277 


[6] 


LSCO, x=012, UD 
LSCO, x=0.26, OVD 


0.8 


1.41 
1.78 


6.84 


290 
463 



Table I: Low-temperature spectral weight in several families of cuprates and several doping (UD=underdopcd, 
OPD=optimally doped, OVD=ovedoped). The data at T — are either extrapolated from the data at T > T c , 
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as in Rcf. [1-3, 6-8], or are taken at T f» T c , as in [4, 5]. Observe that the absolute value of the spectral weight 
determined in Ref. [2, 3] is affected by an error of 15-20 %, which justifies the fact that ilp it is lower in the OPD 
sample than in the UD one. * (The data for YBCO extracted from Ref. [5] refer to the integral of the spectral weight 
up to 0.5 eV. However, the plasma edge reported in this article is still of order of 1 eV.) 



In Table I we reported the values of the spectral weight measured at low temperature as a function of doping. The 
first observation is that the partial spectral weight increases when the doping increases. This behavior is in agreement 
with the general observation that cupratc superconductors are doped Mott insulators [22]. Indeed, these systems 
behave as if the effective number of charge carriers is proportional to the doped holes instead of the total number 
of electrons, so that the spectral weight increases when the system is doped with respect to half filling. It is worth 
noting that this same behavior was observed long ago for the doping dependence of the superfluid density [23] , which 
is also proportional to the partial spectral weight, as we shall discuss below. 

While the doping dependence of the T = spectral weight is in some sense expected for a doped Mott insulator, the 
recent experimental data focused on some new interesting aspects related to the temperature dependence of W(T), 
both in the normal and in the superconducting state. As far as the range of temperature T > T c is concerned, three 
main observations have been made: (i) in all the compounds analyzed in Ref. [1-6] W(ujp, T) shows a T 2 decreasing as 
the temperature increases, but with a coefficient B(ojp) much larger than the one derived in Eq. (2.18) for a standard 
Fermi liquid; (ii) when the cut-off frequency ojm is lower than the plasma edge top the coefficient B(ujm) becomes 
even larger; (iii) in underdoped compounds the opening of the pseudogap does not affect the restricted partial weight, 
in particular it does not lead to any suppression of the spectral weight, as one could expect as a consequence of a 
(pseudo)gap opening. 

This last observation is intimately related to the typical behavior of the spectral weight below the SC transition, 
as described by the BCS relation (2.24). Indeed, if one assumes that the pseudogap opening is due to some form a 
preformed Cooper pairing, then one could also expect to see some decrease of the spectral weight, due to the same 
mechanism which leads to the decrease of W across the BCS transition. The intriguing question arises how to reconcile 
the opening of the the pseudogap with an almost conventional temperature dependence of the spectral weight, where, 
however, the effect of the correlations appears in a renormalization of the coefficient B of Eq. (2.18). To be more 
specific, one can extract from the experiments the coefficient B exp as: 

w exp (n,T) = w exp (n,o) - B exp (n)T 2 . (2.28) 

According to the relations (2.18) and (2.26) for = u>m = ujp one can relate B exp to the value B = w 2 /12D of 
Eq. (2.18) obtained within a tight-binding model: 

ire 2 n 2 E C B 

Bexp ~ ~24D ~ — (2 ' 29) 

We can then extract from B exp an effective hopping parameter tx = tt 2 /48B corresponding to the thermal variation 
of the spectral weight. As one can see in Table II, the tx values extracted in this way are significantly smaller than 
the ones obtained by ARPES measurements of the Fermi surface of cuprates, which are usually of the order of t ~ 300 
meV, leading to a temperature variation B ~ 1/tx larger than expected in a standard Fermi liquid. Even in samples 
were the quadratic temperature dependence has not been clearly observed, the large temperature variation is made 
apparent by analyzing the relative spectral weight variations AW/W = [W(0) — W(300K)]/W(0), reported in the 
last column of Table II. Observe indeed that in the usual (non-interacting) tight-binding model AW/W never exceeds 
the value of about 0.5% [12]. It is worth noting that in the first works of Molegraaf et al. [1] and Santander-Syro et 
al. [2] on the spectral- weight behavior on cuprates the attention was focused mainly on the absence of spectral-weight 
decrease below the pseudogap temperature, but not on the anomalously large values of B measured in these BSCCO 
samples. This issue was addressed for the first time from the theoretical point of view in Rcf. [11]. More recently 
Ortolani et al. [6] clearly stated this problem in their analysis of the optical data of LSCO compounds. As we shall 
discuss in the following Sections, motivated by the observation of the anomalous large value of B we investigated 
indeed the possibility to reconcile the large slope of W(T) with a DDW model for the pseudogap phase of cuprates. 
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Ref. 


Compound 


B (cV- 1 ) 


tr (meV) 


AW/W(0) 


[1] 


Bi-2212, T c =66 K, UD 


9.3 


21 


0.021 




Bi-2212, T c =88 K, OVD 


11.3 


18 


0.022 


[2, 3, 7] 


Bi-2212, T c =70 K, UD 


10.8 


18 


0.042 




Bi-2212, T c =80 K, OPD 


10.8 


18 


0.053 




Bi-2212, T c =63 K, OVD 


16.4 


12 


0.053 


[8] 


Bi-2223, T c =110 K, OPD 


4.2 


47 


0.018 


[6] 


LSCO, x=012 


8.3 


24 


0.042 




LSCO, x=0.26 


10.4 


19 


0.026 



Table II: Slope B of the T 2 temperature variation of the spectral weight in several families of cuprates. We report 
here the value obtained integrating the spectral weight up to the cut-off ojm — uop as reported in Table I for the 
same compounds. The 'thermal' estimate of the hopping tj- is evaluated using a flat band dispersion, see discussion 
below Eq. (2.29). Last column: relative variation of the spectral weight between T = and T = 300 K, AW/W = 
[W(0) - W(300K)]/W(0). 



Finally, we describe the evolution of spectral weight in the superconducting state. When the SC order is established, 
part of the low-frequency spectral weight moves into the SC response, represented by a singular contribution to the 
optical conductivity: 



SO' I \ 

sinqular l^V 



e 2 n s 



n6(w) + - 



(2.30) 



where n s is the superfluid density, related to the penetration depth A by the usual relation 1/A 2 = 4ne 2 n s /m. As a 
consequence, below T c the real part of the conductivity reads: 

2 

Recr(cj) = -6(u)) + a req , (2-31) 

m 

where a reg is the regular part of the conductivity below T c . As far as the sum rule (2.9) is concerned, one sees that the 
condensate contribution is always included in W, and the effect of the SC transition is already taken into account by 
the general formula (2.13) through the change in the occupation number. For example, in the BCS theory the spectral 
weight is expected to decrease below T c , see Eq. (2.24). However, from the experimental point of view the definition 
(2.31) is useful to correctly determine the spectral weight W by taking into account also the condensate contribution. 
Thus, one usually integrates the spectral weight of the regular part up to u sc , W reg (uj sc ) = J + c Reoy eg (f2)cL->, obtained 
from the finite-frequency measurement of the conductivity, and estimate the superfluid weight as W s = ire 2 n s /m, 
usually using the slope of the inverse imaginary part of o-f^ gu i ar , see Eq. (2.30). Then the value W = W S + W,. e9 (w sc ) 
can be compared with the theoretical prediction, which usually concerns the spectral weight integrated up to the 
cut-off Use = cjm- In particular, the BCS model (2.24) predicts a decrease of W as T decrease in the SC state, to 
be contrasted to the increase expected in the normal state according to Eq. (2.18). Even though all the procedure 
described above can be affected by a significant error, due mainly to the exact determination of the superfluid weight, 
nonetheless it was reported by two different group that in an undcrdoped and an optimally-doped Bi-2212 compound 
[1] and in an underdoped Bi-2212 film [2] there is an effective increase of the spectral weight in the SC state with 
respect to the value in the normal state. More specifically, it was shown that by integrating W reg up to the cut-off 
u>M ~wp W(T) is still increasing as T decreases, with a temperature variation even larger than the value B found in 
the normal state, see Eq. (2.28). This result has been recently confirmed in a slightly underdoped Bi-2223 sample [8]. 
Later on, Santandcr et al. showed that in the optimally-doped sample the spectral weight is almost constant in the 
SC state, with W(T < T c ) w W(T C ) [3], as it seems to be the case also for LSCO [6]. Finally, a more conventional 
result has been observed in an overdoped Bi-2212 film [7], showing that the anomalous spectral- weight behavior can 
be eventually doping dependent. However, these results has been questioned in Ref. [5], where the optical data of 
an optimally-doped YBCO sample and slightly-underdoped Bi-2212 sample up to a larger cut-off of ~ 1.5 eV were 
considered. According to the analysis of Boris et al., the spectral weight of YBCO and Bi-2212 in this frequency range 
keeps constant in the normal state and decrease when the SC state is formed. Thus, the issue of the spectral-weight 
behavior below T c remains at the moment quite controversial [9, 10], and seems to be related mainly to the analysis of 
the optical data rather than to an effective discrepancy between them. In particular, it seems to us that in Ref. [5] it 
is mainly questioned the choice of the plasma edge wp ~ 1 cV as the proper cut-off u« for the analysis of optical data. 
However, this is the most plausible choice from the theoretical point of view, because extending the spectral-weight 
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integration to intcrband processes one naturally recovers the full f-sum rule (2.14), which is just a constant and does 
not provide any significant information on the interactions acting on the system. 

III. OPTICAL CONDUCTIVITY SUM RULE FOR A REDUCED d-DENSITY WAVE MODEL 

While it is clear that any successful theory of HTSC should be able to explain the above-mentioned experimental 
facts, despite the 20 year jubilee since the discovery of HTSC, there is no consensus yet neither on the general theory 
of HTSC, nor on the particular explanation of these facts. 

Although our goal is to present here the latest theoretical results for the optical-conductivity sum rule within the 
d-density wave paradigm, we would like to mention some other works related to this issue in the recent literature. A 
first class of papers is related to the spectral-weight below the superconducting critical temperature. For example, 
the possibility of a spectral-weight change below T c in terms of the lowering of the in-plane kinetic energy has been 
analyzed in Refs. [24-27]. In Ref. [24] a model with "occupation modulated" hopping terms was considered. The 
reduction of the kinetic energy at T c was attributed in Refs. [26, 27]) to the transition from a phase-incoherent Cooper 
pair motion in the pseudogap regime above T c (see, e.g. review [28]) to a phase coherent motion at T c , while in Ref. [25] 
a model with a frequency dependent scattering rate was used. More recently, the optical-conductivity sum rule was 
analyzed for a model with electrons in the finite band coupled to a single Einstein oscillator [29] and in the nearly 
antifcrromagnctic Fermi liquid model [30]. 

As far as the anomalous behavior of the spectral weight in the normal-state is concerned, a possible interpretation 
was recently given by Toschi et al. [22]. In this work the authors evaluate the spectral- weight behavior within the 
dynamical mean-field analysis of the repulsive Hubbard model, and they show that the large slope B (sec Eqs. (2.18) 
and (2.28)) is a consequence of the bare-bandwidth renormalization B ~ l/Zt where Z — * as the Mott-insulator 
is approached. At the same time the T = value of W is much less affected by the Z renormalization, due to the 
contribution of the so-called mid-infrared processes, which occur below up and partly compensate for the reduction 
of the quasi-particle weight [22] . 

Finally, the behavior of the spectral weight in the DDW state was addressed specifically in a series of recent papers, 
Ref. [11-13, 31, 32]. While the issue of the vertex corrections, discussed below, is the most important ingredient of 
the Refs. [11-13], in Ref. [31] they were ignored and the main attention was paid to the analysis of the frequency 
dependence of o~(uj) both in the presence and in the absence of the next-nearest-neighbor hopping. As we will see in 
Sec. IV, the results of Ref. [31] are particularly useful for the discussion of interband optical transitions when the vertex 
correction becomes unimportant. At the same time, Ref. [32] investigates the consequences of a momentum-dependent 
scattering time on the optical spectra, neglecting again the contribution of vertex corrections. 



Since 2000 the most commonly used name for the state that we consider in the present article is d-density wave 
(DDW) state [33]. It is considered as one among competing theoretical approaches that are attempting to explain 
the pseudogap behavior in cuprates [21]. However, we do not intend to analyze here to which extent the DDW 
model can explain the vast body of the experimental data accumulated during the last 20 years. Instead of this we 
will concentrate on rather important theoretical subtleties of the DDW model that are often overlooked and their 
implications for the optical conductivity. Yet before writing down the model Hamiltonian we find useful to remind 
the reader of the rather long history of the model. 

The DDW or, as it was called originally, the orbital antiferromagnetic state was introduced in solid state physics 
in Ref. [34] . The discovery of HTSC and searches for a model that would describe the CuC>2 planes without phonons 
lead to the investigation of the Hciscnbcrg-Hubbard model [35-38] and to the rediscovery of this state which at this 
time was called a staggered flux phase. The phase with circulating orbital currents turned out to be an attractive 
idea and there are hundreds of papers (see, for instance, very few of them in Refs. [39-43]) that exploit and develop 
it. These studies showed that the DDW state can be obtained when the interaction term in (2.1), for example, is 
Heisenberg 



A. Unconventional charge density waves 




(3.1) 



<ij> 



or Hubbard-likc 




(3.2) 
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The sum in Eq. (3.1) runs over nearest neighbors sites i and j, a is the vector of Pauli matrices and V(ri — rj) in 
Eq. (3.2) is a generic density-density interaction. For example, the repulsive Hubbard model corresponds to the case 
V{vi — Tj) = U for i = j. 

In the mean-field treatment the expectation value of any two operators, i.e. the order parameter, may be nonzero. 
In particular, if the expectation value of two operators with the same spin projection is nonzero we deal with a 
charge-density wave (CDW) , which can have an internal symmetry if the two operators in the order parameter belong 
to different sites. CDW with an internal symmetry are called unconventional charge-density waves to distinguish 
them from conventional ones, where the electron and hole sit on the same site. An other possibility is that the 
order parameter describes microscopic currents running on few sites of the lattice, as it is the case for the flux phase 
[35-38] or for the current patterns proposed by Varma [39]. It is worth noting that strictly speaking the flux phase 
does not present a charge modulation, even though a gap is present in the excitation spectrum. Indeed, the same 
phcnomcnological spectrum of the DDW state has been proposed as emerging in cuprates due to the tendency of 
the system to form charge order. This attitude was for example considered in Ref. [42], where the DDW state 
emerges as the consequence of the proximity of the system to a quantum critical point [39, 44]. More generally, 
unconventional charge density waves have been observed in many other systems, as for example organic conductors 
[43] or dichalcogenides [45, 46], but we will restrict our discussion here to cuprate superconductors. In the following 
discussion it will be clear that the crucial point in the evaluation of the optical conductivity for a state with DDW 
or, more general, with a CDW with internal symmetry, is the non-trivial dependence of the electronic self-energy on 
the momentum. For this reason the results reported below are generic to many CDW systems, even though we will 
consider specifically the case of DDW order in cuprate superconductors. 



B. Reduced d-density wave model 

In the DDW case, the general interaction term (3.2) is replaced by an approximate reduced interaction Hamiltonian: 

Hl = ~ 277 ^ Wd ( k ) Wd ( k ') C k ( r C k+Q CT c L+Q ( T' C k'<T' , (3.3) 

k.k' 
a,a' 

where Wd (k) = (cos k x a — cos k y a)/2 corresponds to d-wave symmetry with respect to the discrete rotation group of the 
square lattice. By defining iD^ = — (Vq/N) X^ko- w rf(k) < C k+Q<r c k<r > an d adding the non-interacting Hamiltonian 
Hq from (2.1) we arrive at the mean-field DDW Hamiltonian 

H = 53( £k ~ ^) c L c k<r + iDkcl^Ck+Qa], (3-4) 
k,<r 

where Z?k = DoWd^k) is the gap, known as the DDW gap [33], arising from the formation of the state with circulating 
currents below a characteristic temperature Tudw- In writing Eq. (3.4) we considered the case of t' = in the bare 
band dispersion, so that the nesting condition £k+Q = — £k is satisfied, where Q = (7r/a,7r/a) is the characteristic 
wave-vector at which a particle-hole coupling in the DDW state is considered. However, the following results are 
easily extended to the case t' ^ (see, for example, Refs. [12, 31, 47]). 

The notation is then simplified by halving the Brillouin zone and introducing two-component electron operators 
(the DDW equivalent of Nambu spinors [48] widely used in the theory of superconductivity [17]) 

*L=(cL 4+q,.)' (3-5) 
where cl and c^ a are the Fourier transforms of c\ a and c !(T . The Hamiltonian (3.4) written in terms of x becomes 

RBZ 

H = ^2 xia [ £ k°3 - D k a 2 - fj] Xk<r, (3.6) 

k,cr 

where the sum is taken over the reduced (magnetic) Brillouin zone (RBZ). Thus one can easily see that when D$ ^ 
two quasiparticles excitation branches are formed with dispersion £±,k = — M i ^k, where = ye? + D^. As 
a consequence, at finite doping with respect to half-filling the Fermi surface consists of small pockets around the 
(±7r/2, ±7r/2) points and a finite lcading-cdgc shift LE ps Dq — p, at the M points, which mimic in some sense the arc 
of Fermi surface observed by ARPES measurements in the underdoped phase of cuprates [31, 42]. For this reason, 
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the Hamiltonian (3.4) can be considered as an "effective" model for the normal (i.e. non SC) pseudogap phase of 
cuprates, where SC forms within a standard BCS mechanism, since all the anomalies related to the presence of a 
pseudogap are already included in the Hamiltonian (3.4). This point of view was adopted for example in Ref. [42], 
where the phase diagram of Fig. 2 was derived, in excellent agreement with the experimental data on BSCCO. 

400 
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g 250 

g 
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h? 200 
V 

,_ & 150 
100 
50 


0.1 0.14 0.18 0.22 0.1 0.14 0.18 0.22 

S 8 

FIG. 2: Phase diagram of cuprates [taken from Ref. [42]] derived adding the BCS pairing term (3.12) to the Hamil- 
tonian (3.4). Tddw i s the temperatures where Do forms, T* is the temperature where a leading-edge (LE) shift ap- 
pears in the quasiparticle spectrum at the M points, and T c is the SC temperature, at which the SC gap A forms. 




If one consider the Hamiltonian (3.4) as an effective model, it is straightforward to calculate the current-current 
correlation function and the electrical conductivity for it, because it describes non-interacting quasiparticles. Nev- 
ertheless, one should be careful and take into account that the current operator and the diamagnetic tensor for the 
Hamiltonian (3.4) are different from those given by Eqs. (2.10) and (2.11) for the Hamiltonian (2.1). 

This can be understood by deriving the particle current operator compatible with the charge conservation law and 
with the equations of motion for the operators c and c\ [47, 49-51] 

j( q ' *) = N K C Lq/2 CT C k+q/2 CT - ™g ci_ q/2 /k+Q+q/2,] , (3-7) 

where v£. = —dDk/dki is the DDW gap velocity. The first term of the previous expression relates as usual the particle 
current to the band velocity v£, see Eq. (2.12). The second term, which only appears for non-vanishing Dq, takes into 
account the contribution of the orbital currents to the electrical conductivity, and it emerges when the DDW order 
develops. Observe that in the reduced model (3.4) the term proportional to Dq appears as an additional, temperature 
dependent band, which couples k and k + Q electrons, and as a consequence a corresponding term appears in the 
definition of the current. By rewriting the electric current operator (3.7) using the spinors (3.5), one has 



- RBZ 

Ji(0.*) = ^ExL^( k )Xk«r, (3-8) 

where Vi(k) is a generalized velocity which extends the definition (2.12) to the case of Hamiltonian (3.4) 

Vi(k) = vE ( <T 3 +vg ( a 2 . (3.9) 

Note that this result appears to be rather trivial, if one considers the original paper [35], where the staggered flux 
phase was introduced, and the sum of hopping and staggered flux terms was considered as the "effective" hopping 
term. 

The modification induced by the DDW order parameter on the current operator demands a consistent change in the 
diamagnetic tensor with respect to the definition (2.11). As we discussed in Sec. IIA, after the Peierls substitution 
both the current operator and the diamagnetic tensor can be derived from H(A). according to Eq. (2.2). As a 
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consequence, by performing the Peierls substitution in the reduced model (3.4) not only the current operator but also 
the diamagnetic tensor (ja) is modified, containing an extra term for D ^ [11], 



2^ S [ £ k(4aC kCT ) + iD k {cl a c k+Q(T )\ . (3.10) 



When the operator averages are evaluated (see Rcfs. [11, 12] for the details), one finds that the sum rule for the 
reduced model is: 

W DDW (D T\ 1 RBZ 

Eq. (3.11) was derived using the fact that d XiV v£ = 2ta 2 cos k X:V a (and d X:V v^ = ±(Do/2)a 2 cos k x ^ y a), and it reduces 
to Eq. (2.13) for D = 0. Another way to derive the result (3.11) is to evaluate the commutator of the density jo(q, t) 
and the current ji(q, t), given by Eq. (3.7), as explained in Refs. [11, 12]. 

As we mentioned above, one can consider the superconductivity formation on top of the DDW state as arising from 
a conventional BCS type of pairing. More specifically, we can add to the Hamiltonian (3.6) an additional d-wave 
mean-field pairing term 



H p = ]T[A£ C _ ki c kT + h.c], (3.12) 



where the SC <i-wave gap A k is defined after Eq. (2.24). In the SC phase the occupation numbers of Eq. (3.10) get 
modified, and the spectral weight in the DDW+SC state reads 



W DDW (D,A,T) J_ v \ £+ E+ 
Tre 2 a 2 2VN ^ 



RBZ 



- — tanh tanh 

E, 2T E_ 2T 



(3.13) 



where E±^ = y£,± k + A k is the quasiparticle dispersion in the presence of pairing, and the explicit dependence on 
k has been omitted. 



C. Results for the reduced model 



The numerical calculation of the spectral weight defined by Eqs. (3.11) and (3.13) was done in Rcf. [11] and the 
results at doping S = 0.13 are reported in Fig. 3. As one can see, below Tddw the spectral weight W DDW defined in 
Eq. (3.11) is larger than the value obtained in the absence of the DDW order, represented by the dotted-dashed line in 
Fig. 3. When also SC is established, the spectral weight W DDW (D, A) given by Eq. (3.13) is slightly smaller than in 
the DDW state only, the change being proportional to the ratio A/ Do which is quite small at this doping. In the inset 
of Fig. 3 we plot W(T) as a function of T 2 . One can see that below Tddw the T 2 temperature dependence of W(T) 
is still recovered over a wide range of temperature, but with a slope much larger than in the non-DDW case. All these 
results show some interesting resemblance to the experimental data for cupratcs discussed in Sec. II-C. In particular, 
it turns out that the large spectral-weight variation measured by the experiments of Ref. [1-3, 6-8] between T = 
and T = 300 0. It can be ascribed to existence of DDW order in the normal state. Moreover, when SC is formed 
on top of this DDW state the spectral weight remains almost constant, as observed in Ref. [3]. However, one should 
notice that on the overdoped side, the temperature Tddw of DDW formation decreases (see Fig. 2) so that Tddw 
can be smaller than room temperature, and a change of slope of the spectral weight at Tddw should be observed. 
Up to now, no signature of this effect has been observed in the experiments. 

To develop a better insight into the origin of the spectral-weight increase in the DDW state we calculated in Ref. [12] 
the optical conductivity corresponding to the sum rule (3.11). This can be done using the general definition (2.6) 
where, however, the correlation function (2.5) is modified according to the new current definition (3.7). Let us first 
introduce the Green's function (GF) corresponding to the effective Hamiltonian (3.6): 

G _1 (k, iuj n ) = (ioj n + (i)a - e k cr 3 + D k <7 2 . (3-14) 

Using this GF one can write down an exact expression for the correlation function 

j, RBZ 

Hy(q = 0,iSlm) = -2— Tr[G(k,wj„ + ffi m ) 7i (k)G(k, iw n )7 j (k)], (3.15) 
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FIG. 3: Spectral weight in units of e 2 na 2 /V in the normal state Eq. (2.17) (dot-dashed line), in the DDW state Eq. (3.11) (dashed 
line) and in DDW+SC state Eq. (3.13) (full line). The values of parameters are the same used in Fig. 2 for the doping 8 = 0.13 
(D (T = 0) = 0.92t, A(T = 0) = 0.0644, see Ref. [11, 42] for further details). The critical temperature is marked by the arrow. 
Observe that the decrease of W DDW (D, A, T) below T c is almost negligible. Inset: spectral weight plotted as function of (T/t) 2 . 



where 

7i (k) = ViQt), (3.16) 

with Vi(k) defined in Eq. (3.9) is the bare vertex for our noninteracting theory [52]. One can also convince oneself 
that (3.16) is indeed the appropriate vertex for the reduced Hamiltonian (3.4) by looking at the Ward identity (see 
Eq. (4.4) below) 

which is just the charge-conservation law expressed in terms of the Green's and vertex functions. Then one can state 
that the term ~ giv® in the vertex (3.16) is related to the k-dependence of the DDW gap in the GF (3.14). Using 
the definition (3.15) of the correlation function and adding an isotropic impurity scattering we evaluated in Ref. [12] 
the optical conductivity in the DDW state, see the inset of Fig. 4. Here we report the T = value of the optical 
conductivity in a system without DDW and with DDW for 6 = 0.13, using the same parameter values as in Fig. 3. 
As usual, when DDW order forms the optical conductivity presents both a Drude peak and an absorption at finite 
frequency ~ 2[i due to the optical processes between the two quasiparticle branches £± (see also Ref. [31] for the case 
with t' 7^ 0). Although the Drude peak is reduced by the DDW formation, these interband finite- frequency processes 
lead to an overall increase of the spectral weight with respect to the case without DDW, as one can see in the main 
panel of Fig. 4 where we report N(u) = 2 J? a(uj')duj' . The crossing of N DDW with respect to N at w ~ 2t is again 
an effect of the additional term proportional to vd in the current (3.7) and in the vertex (3.16), which extend above 
the Drude peak the effect of the DDW formation. Indeed, as discussed in detail in Ref. [12], and consistent with the 
results of Valenzuela et al. in Ref. [31], when the contribution proportional to is omitted one finds instead in the 
DDW state a spectral weight lower than in the normal state (in analogy with the behavior across the SC transition 
described by ordinary BCS theory). 



IV. OPTICAL CONDUCTIVITY AND SUM RULE WITHIN A FULL MODEL WITH DDW 

INSTABILITY 

A. Limitations of the reduced model and need of a more microscopic approach 

There is no doubt that the results presented above for the electric current operator (3.7), the optical conductivity 
shown in Fig. 4 and the sum rule (3.11) (and its generalization (3.13) valid in the superconducting state that forms 
on top of the DDW state) are consistent with each other, as far as they are considered within the framework of the 
reduced model (3.4). 
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t-0.25 eV, 5=0.13, T=0 




FIG. 4: Integrated spectral weight in units of e 2 ira? /V for the optical conductivity at T = of a system with 
and without DDW. The values of parameters are the same as in Fig. 3. At low cut-off energy the spectral weight 
in the DDW state is lower than in the normal state, but when intraband excitations are taken into account the spec- 
tral weight of the DDW state exceeds the value of the normal state. Inset: optical conductivity in the two cases. 



One could, however, still question these results from the more general point of view of the validity of the Hamiltonian 
(3.4), viz. one may question whether the low-energy Hamiltonian (3.4) is suitable for the description of the intraband 
excitations that cause the increase of W DDW (D, T) as T decreases. The justification of the Hamiltonian (3.4) can 
only be made using a more microscopical Hamiltonian (2.1). In many ways this resembles the history of the reduced 
BCS Hamiltonian, since the theory of superconductivity was completed only when a more generic Hamiltonian was 
studied (see e.g. Refs. [17, 48]). These analogies between DDW and BCS models were exploited in Ref. [12]. One of 
the main conclusions of Ref. [12] is that the relation between the optical conductivity a{u>) and the sum rule for the full 
model (2.1) can be established by looking at the problem of the vertex corrections in the current-current correlation 
functions. This vertex is determined by solving an integral equation, whose solution is a formidable problem. Its 
zero-frequency solution provided by the Ward identity was used in Ref. [12]. More recently Aristov and Zeyhcr [13] 
extended this analysis to finite frequency and investigated the behavior of sum rule in the full DDW model. 



B. Electrical conductivity and equation for vertex 

The current-current correlation function defined by Eq. (2.5) is expressed in terms of the DDW Green's function 
(3.14), the bare vertex yf and the full vertex function Tj as follows [17]: 

rp FIB Z 

n„( q = o,io m ) = -2— Tr[G(k, jw n ) 7 ?(k)G(k, w„ + in m )r j -(k,*n m )]. (4.1) 

where 7 and V represent the bare and full vertex function, respectively. A careful reader may notice that the vertex 7? 
for a model based on the full Hamiltonian (2.1) is nothing else then Eq. (2.12). Using the matrix notation introduced 
in Sec. Ill it reads 

7 °(k) = (^a 3! a ). (4.2) 

Comparing the bare vertex 7? with the vertex 7$ [see Eqs. (3.16) and (3.9)] of the reduced model, it becomes clear 
that the bare bubble approximation when the full vertex Tj(\s.,u>) is replaced at will by the bare vertex 7j°(k) would 
miss completely the physics of the problem. Thus one has to find the full vertex by solving the integral equation for 
it [12, 13] 

r«(k, in m ) = 7 °(k) [ dpV{p - k)G(p, ixi)Ti(p, iCl m )G(p, ixt + iO ro ), (4.3) 

xi J 

where the current vertex 1^ is a 2 x 2 matrix and the integration over p with the measure dp = dp x dp y / (2ir) 2 runs 
over the full (chemical) Brillouin zone, while q lies within the reduced zone. Here V(p — q) represents the Fourier 
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transform of the potential V(r, — rj) of Eq. (3.2). In Ref. [13] the interaction term of the Hamiltonian Hi nt is 
chosen to be the Heisenberg interaction (3.1) between nearest neighbors, so that in the momentum representation 
V(k) = J(k) = 2J(cosk x a + cosk y a). After analytical continuation iCl m — > to + iO is made, one can consider the 
static limit u> = of Eq. (4.3) to illustrate that if a vertex function I\ satisfies the Ward identity 

W)=-^ (4 . 4) 

it is also a solution of the equation (4.3) for vertex. The crucial assumption in this proof (see Refs. [12, 17]) is that 
the potential V is nonseparable. viz. it depends on the difference p — k. The restrictions on the form of the potential 
V(p, k) necessary to satisfy the Ward identity (4.4) and, respectively, to obey the charge conservation law deserve 
special study (see e.g. Refs. [53, 54]). It was observed in Ref. [13] that an interesting feature of Eq. (4.3) is that, for 
the specific potential J(p — k) considered there, only the odd kernels (like smp x asmk x a) of its decomposition into 
the sum of separable kernels can contribute in the integral over p, while these terms were absent in the reduced model 
(3.3). Finally we note that the relation (4.4) was used above in Eq. (3.17), since for the reduced non-interacting 
model (3.4) it allows to obtain the exact vertex valid for all frequencies [12]. 

The equation for the vertex (4.3) with the potential J(k — p) written above is solved at finite frequency in Ref. [13]. 
The corresponding optical conductivity (2.8) calculated on the base of the full bubble (4.1) reads 

<t(w) = -Im fdkf dx x dx 2 np{Xl) ~ nF{X2 l Tr[ 7t (k)A(k, xJE^k, iu)A(k, x 2 )\, (4.5) 

where the integration is over RBZ and 2wiA(u),k) = G^(w,k) — Ga(uj,k) is the spectral function expressed via the 
difference of advanced and retarded GF. Note that to regularize the electrical conductivity in the limit to — > one has 
to consider the scattering by impurities which is included in the simplest form by introducing a finite quasiparticlc 
lifetime, r qp . The expression for <r(u>) becomes rather complicated when the solution of (4.3) is substituted there, but 
in the limits u) — ► and LOT qp 3> 1 it can be significantly simplified. These two limits are considered below. 



C. Conductivity for c 

According to the results of Refs. [12, 13], in the limit u — > the conductivity <t(uj) reduces to the "naive" expression 

a(cj) naive = -Im fdk f dx x dx 2 np{Xl) ~ np{X2) Tr[r a (k, 0)A(k, x^k, 0)A(k, x 2 )], (4.6) 
to J J u> — Xi + x 2 + «lJ 

with the static vertex r^k, 0). The dc limit of o~{oj) na i ve was obtained also in Ref. [51], where the link between the 
present problem and the approach developed by Langer [55] was pointed out. 

Comparing the LHS of Eqs. (4.4) and (3.17) which enter, respectively, in the current-current correlation function 
for DDW model (3.15) and in Eq. (4.6), one can notice that o~{oj) na ive coincides with the exact <j(uj) derived from 
Eq. (3.15) for the reduced DDW model. This coincidence was the main motivation for the work of Ref. [12], where 
the conductivity a(u)) na i ve , which is just an approximation for the full model (2.1) with H int given by Eqs. (3.1), was 
extrapolated to higher tu and studied in detail. 

The analysis of Ref. [13] clarifies that Eq. (4.6) does not remain valid for higher co. Considering the general shape 
of the optical conductivity depicted in Fig. 4, one finds that: (i) as far as the Drude peak is concerned, for tuT qp 3> 1 
the current vertex r^(k, u>) changes from its to = value Eq. (4.4) to its bare value Eq. (4.2). Accordingly, the value of 
<j{u>) already near the tail of the Drude peak may well be estimated using Eq. (4.5) with two bare vertices 7°(k); (ii) 
as far as intraband transitions are concerned, the vertex correction enhance cr(uj) only near the edge of the intraband 
optical absorption. Thus the expression (4.6) overestimates the value of o~{u>) both in the Drude part, for u>r qp 1, 
and in the region of intraband transitions. A somewhat counterintuitive result of Ref. [13] is that the crossover from 
the full vertex Ti(k,u>) to the bare vertex 7?(k) occurs at u ~ r qp rather than at the energy scale Do related to the 
formation of the DDW order. It would be interesting to study how this result depends on the choice of the microscopic 
Hint ■ We note also the influence of the full vertex Ti (k, uj) discussed here is obtained under the assumption that the 
impurity scattering rate \/2r qp is isotropic. As shown recently in Ref. [32], these results may change drastically if 
l/2r qp acquires a k-dependence. 



D. Behavior of the spectral weight 



The temperature dependence of the spectral weight W(T) obtained in Ref. [13] turns out to be opposite to the 
behavior of W(D,T) considered in Sec. Ill, viz. in the DDW state it decreases as the temperature T decreases and 
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the DDW gap D(T) opens. This decrease follows from the fact that for Do > r^ 1 the vertex r^k, lu) ~ 7°(k), so 
that the weight of intraband excitations is smaller than in the reduced model, where it is crucial for producing an 
increase of W DDW (D, T) that the vertex 7j(k) [56] is effective for all frequencies u> (see Sec. III). Observe that this 
result should be expected according to the fact that the kinetic energy increases when the DDW state is formed. As 
was explained in detail in Ref. [12], when vertex corrections are properly included in the current-current correlation 
function the resulting sum rule is given by Eq. (2.13) with a modified occupation number due to the DDW transition, 
so that W(T) decreases, in analogy to what W(A,T) given by Eq. (2.24) does in the superconducting state. 

The results of Ref. [13] also show that in the presence of the full vertex the absolute value of the spectral weight 
significantly increases with respect to the estimate made with the bare vertex 7 l °(k). Moreover, this increase is seen 
already in the normal state above T DDW , so that one can consider that even when W(T) starts to decrease below 
Tddw, the absolute value of W(T) is still well above the value of the weight estimated with the bare vertex 7°(k). 
This result should be contrasted with the analysis of Ref. [22, 30], where instead it was shown that the interaction 
decreases the absolute value of the normal-state spectral weight with respect to the non-interacting case, mainly as 
a consequence of the quasiparticle renormalization which is instead absent Ref. [13]. Thus, further investigation is 
required to understand whether this discrepancy is related simply to the choice of different microscopic models or it 
is caused by the approximations used to evaluate the spectral-weight behavior. 

V. OPEN QUESTIONS AND RELATIONSHIP TO OTHER SYSTEMS 

Our purpose here was not to convince our reader that the DDW state is the best candidate to explain the pseudogap 
behavior observed in HTSC. Neither the reduced nor the full model show full agreement with existing experimental 
optical-conductivity data, due in particular to the lack of a satisfactory description of the shape of the optical spectra 
observed experimentally. In particular, no clear signature of the interband processes characteristic of the DDW state 
has been observed in cuprates [20] , even though some signatures of finite-frequency peaks have been reported recently 
[57], which can still be attributed to some kind of charge ordering [58]. Moreover, while the reduced DDW model 
shows at least an interesting feature that agrees with the experiments, i.e. an increasing spectral weight in the non 
SC state, an attempt to generalize it and consider a more microscopic model destroys this behavior. Certainly more 
work has to be done to understand the link between reduced and full models. Nevertheless we hope that the results 
discussed here for the restricted optical-conductivity sum rule in the DDW state may be useful for the investigation 
of this sum rule in more sophisticated models leading to unconventional charge-density wave ordering or circulating 
currents (see e.g. [39]). Since the self-energy in these models is also k-dependent, to study the transport properties 
consistently one should definitely include vertex corrections. Finally we mention that the results presented here may 
be useful not only to cuprates, where they can be only partly applied, but also to other materials displaying a true 
fc-space modulated CDW, as for example organic conductors [43] and dichalcogenide materials [45, 46]. 
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